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Abstract 

The hghter superpartner of the top quark {ti) may be considerably hghter 
than squarks of the first two generations, and hence may be accessible to 
Tevatron collider searches, even if the other squarks and the gluino are too 
heavy. For the range of m^^ of interest at the Tevatron, the ti decays to a 
chargino via ti — > hWi if this is kinematically allowed; otherwise, the flavor- 
changing loop decay ti — > cZi dominates. In the latter case, ti production is 
signalled by jet(s) plus I^j, events. If, instead, the chargino decay is allowed, 
1£ + b-jets + I^rp and + jet{s) + events from leptonic decays of Wi 

provide the most promising signals. We perform detailed simulations for each 
of these signals using ISAJET 7.07 and devise cuts to enable the extraction of 
each of these signals above standard model backgrounds from vector boson or 
top quark production. With an integrated luminosity of 100 pb~^, experiments 
should be able to probe top squark masses up to ~ 100 GeV in each of these 
channels; the detection of the signal in the single lepton channel, however, 
requires reasonable capability to tag displaced b decay vertices in the central 
region. 
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I. INTRODUCTION 



The observation that precision measurements of gauge couphngs at scale Mz by the LEP 
experiments are consistent with the simplest supersymmetric ^ SU(5) grand unification 
(but not with minimal non-supersymmetric SU(5)) has led many authors to seriously re- 
examine the expectations for sparticle masses 0| within this framework. More recently, 
several studies of the collider phenomenology of supergravity models have also appeared . 
On the experimental side, direct constraints on such models come from non-observation of 
super-partners at colliding beam experiments. For instance, the non-observation of multi- 
jet plus missing transverse momentum {I^j) events above expected standard model (SM) 
background levels has led the CDF and DO collaborations to conclude 

m-g, rrig > 100 - 150 GeV, (1) 

where the considerable range in the bound is due to the dependence of the missing cross 
section on the parameters [0 of the Minimal Supersymmetric Model (MSSM). In obtaining 
these mass limits, it is usually assumed that there exist ten or twelve types of mass degenerate 
squarks. This is usually justified by appealing to the framework of minimal supergravity 
models, where the sfermions are all expected to be degenerate at an ultra-high unification 
scale; this degeneracy is broken when these MS mass parameters are evolved down to the 
electroweak scale relevant for phenomenological analyses. For the squarks of the first two 
generations, this evolution is dominantly governed by their (common) QCD interactions so 
that their masses are split only by the relatively small electroweak interactions. As a result, 
these remain essentially degenerate with a mass rrig. 

In contrast, the masses of third generation squarks obtain substantial contributions 
also from the large top family Yukawa interactions, which also induce a substantial mixing 
between the and tR squarks. These interactions reduce the diagonal masses relative 
to ruq in much the same way that they drive the Higgs mass squared to negative values, 
resulting in the breaking of electroweak symmetry. At the weak scale, the soft-breaking 



masses for t-squarks can be written in the form [10 



^Ir = ^1 ~ ^ft^l^t + ml (2) 

Here, ft = ^ v^"^^ + "^^^ ifb = ^/ V^'^ + V^) are top (bottom) family Yukawa cou- 
plings, mo is the universal sfermion mass at the unification scale and f, are dimensionless 
parameters ~ 0.1 that can be numerically computed within a model. The term involving the 
bottom Yukawa coupling is negligible unless tan/3(= ^) is close to mt/rrih [|ll|]- Assuming 
squarks are heavier than top quarks, we expect the mass ordering m^^ < mi^ < nig. As 
already noted, tL and tR are mixed by their Yukawa interactions. The t-squark mass squared 
matrix takes the form, 



m?^ + mf - M| cos 2/3(-i + | sin^ Ow) ~mt{At - /i cot (3) 

-mt{At- fx cot P) mf + mf + cos2P{^sm^ 9w) 



(3) 

where At is the trilinear soft SUSY breaking scalar coupling evaluated at the electroweak 
scale and fi is the superpotential Higgs mixing term. The mass eigenstates [ii and ^2) can 
be readily obtained by diagonalizing this matrix. For the eigenvalues, we have 
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<,2 = ^("^i + "^D + cos 2(3 + ml 



iK-my + ^(8M^-5M|) 



+ m^(A -/icot/3)H . (4) 



The corresponding mass eigenstates are 



ti = cos 6 ft L — sin Oft ji, 

t2 = sin 6tiL + cos OttR, (5) 



where the mixing angle 6t is given by 

m? + - M|cos2/3(-i + |sin2 6'vy) - m? 

tan 0t = ' ^- ^-^ ^ . (6) 

-mt[At — /icot p) 



The freedom to adjust the soft-SUSY breaking parameter At enables us to fix the lighter 
top squark (the stop, or ti) mass essentially independently of the masses of other squarks. 
We see that ti, which has a mass smaller than even m^^, can indeed be much lighter than 
the first two generations of squarks. In fact, ti can easily be as light as ~ 50 GeV even 
if other squarks and gluinos have masses of several hundred GeV. Although this is not 
directly relevant to our discussion, we also mention that SU(2) gauge invariance requires 
that m^^ = mi^ (aside from D-terms) whereas m^^ ~ m^ unless tan/3 is very large. For 

moderate values of tan/5, /i and A;,, hL-hn mixing is expected to be small so that and hR 
closely approximate the mass eigenstates. 

Having convinced ourselves that ii can indeed be much lighter than other squarks, it is 
reasonable to ask what current experiments tell us about its mass. The measurement of the 
total width of the Z boson as well as direct searches for squarks at LEP [0 imply |jl3[ that 



mi^ > 45 GeV. (7) 

Top squarks can also be pair produced at pp colliders, via qq and gg fusion diagrams. 
Unlike squarks of the first two generations, there exists no significant cross section contri- 
bution from gluino exchange diagrams, so that a{titi) is determined completely in terms of 
mj^. For stops with a mass in the 50-125 GeV range — accessible to Tevatron experiments — 
the dominant decay mode is expected to be the two body ti bWi, if it is kinematically 
allowed. If instead m^^ < m^ + ttl^^ , then decay via flavor changing loop diagrams is ex- 
pected to dominate 0], in which case ti — cZi. In the latter case stop production, like 
squark production (but without cascade decays), would be signalled only by events with 
jets plus I^rp. When the chargino decay mode is accessible, stop signatures are topologically 
the same as those of t-quarks. We have checked at least within the MSSM framework that 
if m^^ ^ 150 GeV, the decay ti — > bWZi is kinematically allowed only when the two-body 
decay ti bWi is also allowed. Thus the three body decay to real W bosons is never 
phenomenologically relevant for the discussion of t-squark signals at the Tevatron. 



The production of ti-squark pairs at the Tevatron was studied in Ref. |T^, where it 
was shown that o"(titi) ~ (^ — ^)cr(tt), for equal top and stop masses. These authors. 
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using parton level Monte Carlo programs, attempted to translate the CDF ^ squark bound 
Eq.(l) to a limit on the mass of ii, under the assumption that ii cZi. They concluded 
that the CDF data imply a stronger bound than the LEP limit in Eq. (4), but only if the 
lightest neutralino, also assumed to be the lightest SUSY particle (LSP), has a mass of just 
a few GeV. For m^^ > 10 GeV, their analysis showed that it is entirely possible that ti 
could have escaped detection if its mass was just above the LEP bound. These authors 
also studied the dilepton signal from the decay ii bWi, Wi — * iuZi. Since their purpose 
was to study whether it would be identifiable in the CDF search for top quarks, they used 
the corresponding CDF cuts. They found that the stop signal would be dwarfed by the 
corresponding signal from t-quark production. 

Since that time, the CDF experiment has more than quadrupled its data sample. Fur- 
thermore, the DO experiment has started operation and accumulated a data sample of com- 
parable magnitude. By the end of the current run (Run I B), the combined integrated 
luminosity is expected to be > 100 pb~^; i.e., ^ 25 times larger than that used in the earlier 
analysis [|15|. On the theoretical front, the production and decay patterns of all sparticles 



as given by the MSSM have been incorporated into the event simulation program ISAJET 



||16|| . This should allow for considerably better simulation of the signals than in Ref. [|15|. 

The purpose of this paper is to re-examine the prospects for discovering the lighter 
top squark at the Tevatron in the case where other squarks and gluinos are too heavy to be 
kinematically accessible, regardless of whether the decay ti — > bWi is allowed. In light of the 
fact that the Tevatron experiments will soon have a very large event sample, and considering 
the improvements in theoretical technology since the last study, we felt that a reassessment 
was warranted. For the signal we improve on our earlier study in that we use a more 
sophisticated simulation. We will see later that the inclusion of QCD radiation in ISAJET 
significantly affects our conclusion about the mass reach of the Tevatron, particularly when 
Zi is relatively heavy. In the case when the chargino decay of the stop is kinematically 
allowed, we adopt a different philosophy from our earlier study. Instead of focussing on 
what the Tevatron top quark analysis can tell us about top squarks, we devise cuts to 
separate the SUSY signal from the top signal as well as from other SM backgrounds. We 
then find that, with an integrated luminosity of ~ 100 pb~^, experiments at the Tevatron 
should be able to search for ti as heavy as about 100 GeV, even in this case. 

The rest of this paper is organized as follows. We briefly describe our simulation of stop 
events at the Tevatron in Sec. II. In Sec. Ill, we investigate Tevatron signatures from ti — >■ cZi 
decay, and show regions of the m^^ vs. m^^ parameter space that ought to be accessible 
to Tevatron collider experiments. In Sec. IV, we examine both the single lepton plus jets 
signature as well as the dilepton plus jets signature, when ti — >• bWi is the dominant stop 
decay. Again, we show regions of parameter space explorable via these modes. We conclude 
in Sec. V with a summary of our results along with some general remarks. 



II. TOP SQUARK SIMULATION AT THE TEVATRON 

We use the program ISAJET version 7.07 |l^ to simulate events from the production 
of ti pairs at the Tevatron. Once produced, the stop rapidly decays via ti bWi if this 
decay is kinematically allowed; if not, we assume that it decays via cZi with a branching 
fraction of 100%. In other words, the decay patterns of the stop are determined by sparticle 
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masses and do not depend on the stop or gaugino mixing angles. Furthermore, in this latter 
case, the l^rp signal from stop production is completely determined in terms of m^^ and 
and is independent of any other parameters. If, on the other hand, the stop is heavy 
enough to decay via the chargino mode, its signals obviously depend on the decay patterns 
of the chargino. For any set of input parameters, ISAJET then decays the daughter chargino 
with branching fractions as given by the MSSM. Over a large range of model parameters, 
the IV*-mediated amplitude dominates the decays of the chargino, so that the branching 
fraction for the decay Wi \ivZx is about 11%. In our computations, we will use the 
supergravity motivated choice /i = — as our default value and fix tan/3 = 2. This then 
implies '^fy^ — 2m^^ \'^9^ which together with the stop mass fixes the kinematics of 
the events. Effects from radiation off initial and final state partons, fragmentation of the 
c or 6 daughters of the stop, final state hadronization and underlying event activity are 
incorporated in ISAJET. 

We have modelled the experimental conditions at the Tevatron by incorporating a toy 
calorimeter with segmentation Ar] x = 0.1 x 0.087 and extending to \!]\ = 4 in our 
simulation. We have assumed an energy resolution of 50%/ \/Et {15%/\/Et) for the hadronic 
(electromagnetic) calorimeter. Jets are defined to be hadron clusters with > 15 GeV 
within a cone of AR = v^Arp~+~A0^ = 0.7 and |?7j| < 3.5. We consider an electron (muon) 
to be isolated if pri^) > 8 GeV {pt{ij) > 5 GeV) and the hadronic energy in a cone with 
AR = 0.4 about the lepton does not exceed mm(^£'r(^),4 GeV). Non-isolated electrons are 
included as part of the accompanying hadron cluster. 



III. SEARCH FOR STOP DECAY TO CHARM PLUS NEUTRALINO 

When the chargino is heavier than (m^^ — rrih), the ti in the mass range of interest 
essentially decays via ti cZ\ so that stop pair production would be signalled by 
events with charm jet(s) together with In our analysis of this signature, we have 
imposed the following CDF-inspired cuts on the signal: 

• (z) We require at least two jets in each event with at least one jet in the central region, 
\r]j\ < 1. All jets are required to be separated by at least 30° in azimuth from J^rp. 

• (a) If Tij = 2, we further require A0(ji,j2) < 150°. 

• {in) We veto events containing leptons (from the c-jet) with pt{1)> 10 GeV to reduce 
the background from W ^ iu {i = e or fi). We have checked that this leads to very 
little loss of signal. 

• (iv) We require I^j. > 50 GeV to reduce backgrounds from QCD heavy flavours and 
mismeasured jets. 

SM backgrounds from heavy flavor production (cc and bb) and multi-jet production 
(with substantial jet energy mismeasurement) are very dependent on a detailed detector 
simulation, but are expected to be small, given a large enough I^j, cut and lepton veto. 
The major backgrounds we consider here are from vector boson production, and include (a) 
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Z vv production, (b) W tu production (where the hadronically decaying r can be one 
of the jets), and (c) W iu, where extra jets come from initial state QCD radiation. 

The I^j, distribution, before any cuts, from stop events at app coUider with y/s = 1.8 TeV 
is illustrated in Fig. 1 for two representative choices of ii and Zi masses: (a) m^-^ = 85 GeV, 
— 20 GeV, and (6) m^^ = 125 GeV, = 40 GeV. Also shown are the corresponding 



m 



distribution from the W and Z backgrounds listed above which have been computed using 



ISAJET. We have used the Set I structure function of Eichten et. al. in our computation. 
As expected, the W ^ background peaks at ~ Mt^/2, while the spectra from the 
other backgrounds are softer. We see that the J^j, > 50 GeV cut (iv) is very effective in 
cutting these backgrounds with a factor of ~ 2 loss of signal. Nevertheless, even after cuts 
[i-iv), the W —>■ T {Z —>■ vv) background is 23 pb (12 pb) in contrast to the signal cross 
section of 9.2 ph and 1.9 pb for our representative cases introduced above. The lepton veto 
reduces the W ^ background to negligible levels. We have also checked that Z rf is 
small compared to the other backgrounds. 

To facilitate further separation of the stop signal from SM background, we have studied 
the correlation between A0, the transverse plane opening angle between I^j, and the nearest 
jet, and pT(fast jet) for the I^j, event sample. We have illustrated this by scatter plots in 
Fig. 2 for the two signal cases and for the dominant backgrounds. We see that by further 
requiring, 

• {v) Prifast jet) > 80 GeV for A(j) > 90°; else prifast jet) > 50 GeV, 

we are able to substantially reduce the main SM backgrounds, with relatively modest loss of 
signal (particularly for the ti(125 GeV) case, where the signal was indeed small compared 
to the backgrounds). 

The resulting cross section from titi production at the Tevatron, after cuts (i-v), is 
illustrated by the contour plot in the m^^ vs. m-^^ plane in Fig. 3. Also shown are the 
SM backgrounds after these same cuts. We note that although for certain regions of this 
parameter plane the chargino decay mode of the stop will necessarily be kinematically open 
within the MSSM framework, for the purposes of this figure, we have assumed that the 
charginos are too heavy to be produced via decays of the stops, and used the FORCE 
command in ISAJET to decay the stop via ti cZ\. We terminate the contours for stop 
masses where the three body decay t\ — >■ bWZi becomes accessible, since then, the branching 
fraction for the two body mode cZi mode is no longer unity. Several features of Fig. 3 are 
worth noting. 



1. We see that the signal cross section after all cuts exceeds 1 pb for m^^ as large as 
~ 100 — 120 GeV even if Zi is as heavy as 50 GeV. For this cross section contour, as 
many as 40 stop events may already be present in the accumulated data sample of the 
Tevatron. 

2. The dominant SM background source is W —>■ tu production, where the hadronically 
decaying r is counted as one of the jets. Since r jets almost always have a charged 
multiplicity of 1 or 3, it should be possible to discriminate the bulk of these events 
from the signal by vetoing low charged multiplicity jets. We show this background 
to allow the reader to assess the r-jet discrimination factor which is necessary for 
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sufficient background rejection. If a rejection 80% ofW^r events is achieved without 
substantial loss of signal, the stop signal exceeds combined background for m^^ ^100 
GeV, and m^^ < 30 - 40 GeV. 

3. We note that the Z ^ vv background can be reliably subtracted if a large enough data 
sample is accumulated for directly measuring high Z bosons decaying to leptons, 
and using the Z branching ratios measured at LEP. 



4. Unlike as in Ref. [0, where it was concluded that the stop signal was very small for 
m^^ > 20 GeV, we see that the signal is quite robust even if Zx is heavy. In fact, 

we have checked that there is a significant signal level even if Z^ is nearly degenerate 
with ti — for instance, for m^^ = m^^ + 5 GeV = 105 GeV, the cross section after all 
cuts is 0.25 p6, which is more than 2% of the total stop production cross section of 
11 "ph. The sharp fall-off of the cross section with increasing m^^ in the earlier parton 
level simulation occured because the daughter charm quarks rapidly became too soft 
to satisfy the jet requirements. In the more realistic ISAJET simulation, the stop 
pair can be produced with a substantial p-r since it can recoil against jets from QCD 
radiation. For m^^ ~ '^Zi' stops give most of their momentum to the heavy Zi, 

so that I^j, ^ prititi) and the jets come from QCD radiation as, for instance, in the 
Z —>■ uj? background. 

5. If, as anticipated, the Tevatron experiments accumulate around 100 pb~^ of integrated 
luminosity by the current run, we may anticipate at least 50 or more signal events in 
the data for m^^ ^ 100 GeV even if Zi is relatively heavy. Although the background 
level, as read off from Fig. 3, would be around 200 events (assuming aW ^ t rejection 
of 80%), we should keep in mind that the actual background, after Z ^ vv events are 
subtracted, will be considerably smaller. A several standard deviation signal should 
be possible for stops as heavy as 100-125 GeV. A precise evaluation of the Tevatron 
reach, which should take into account detector-dependent backgrounds, is beyond the 
scope of this study. 

Up to this point, we have made no use of the fact that the signal always contains c-quark 
jets. It is clear that SM backgrounds would be considerably reduced if it were to be possible 
to tag at least one of the c-quarks. This led us to consider the possibility of using a muon 
from the semi-leptonic decay of one of the c quarks as a tag. The signal would then consist 
of /i -|- (rij > 2) -|- [J^j. > 50 GeV), where the muon is within a cone of AR = 0.4 about one 
of the jets. For these muon-tagged events we require, in addition to {i-iv) above, that 

• PT{fJ')> 3 GeV for the muon to be identifiable, and 

• either A(f){jnear,$T)< 90° or pT{3fast)> 50 GeV. 

The signal cross section contours, with these cuts, are shown in Fig. 4 together with 
estimates of backgrounds from W —>■ tu fiuui', W ^ fiu and Z ^ uu + cc or bb processes. 
To estimate these, we have generated 14K (130K) W (Z) events of each type, and find that 
just 8, 5 and 6 events, respectively, pass our cuts. It should also be kept in mind that 
ISAJET does not include the full 2 — > 3 matrix elements for Zee or Zbb production; in our 
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simulation, these events come from radiation of initial state gluons followed by splitting into 
bb or cc pairs. 

We see from Fig. 4 that even with an integrated luminosity of 100 pb~^, there will be 
only 5-10 tagged events for sparticle masses such that the signal shown in Fig. 3 might 
be difficult to observe above the background. We further see that the 0.2 pb contour in Fig. 
4 (which roughly maps out the region where S/B > 1) essentially tracks the contour where 
the I^rp signal in Fig. 3 equals the Z — > z/z/ background. We thus conclude that the use of 
muon tagging does not extend the reach of the Tevatron to discover stop squarks in the 
channel. Nevertheless, observation of /x-tagged I^j, events could be important, even at low 
rates: non-observation of these events at expected rates would rule out a fi — cZi decay 
hypothesis. 

IV. SEARCH FOR STOPS DECAYING TO CHARGINOS PLUS BOTTOM 

If m^^ > -|- rrii,, the stop decays via the two-body chargino mode with a branching 
fraction of nearly 100%. Charginos with ^ 100 GeV typically decay via the three-body 

mode Wi f f'Zi (/ = g or i) mediated by a virtual W, qi or II. Since sfermion masses are 
likely to be larger than M^y, the W* mediated decays dominate over large ranges of SUSY 
parameters, so that the branching fractions for chargino decays are, generally speaking, 
close to those of the W-hoson. It is only when the lightest neutralino is dominantly a U{1) 
gaugino that the WWiZi coupling is dynamically suppressed; then the sfermion- mediated 
amplitudes become especially important. If, in addition, sleptons are also considerably 
lighter than squarks as is the case in the no scale limit of supergravity models, the branching 
fractions for chargino leptonic decays can be significantly enhanced. Here, we have mainly 
confined our attention to the more typical case where iy*-mediated decays indeed dominate. 
It is, however, worth remembering that the leptonic signals discussed here may be enhanced 
in some regions of parameter space. 

If both the charginos in a stop pair production event decay hadronically, the signal is 
qualitatively the same as that in Sec. Ill: i.e., consisting of n — jet + I^j, events. Since, in 
the present case, the Zi is always produced via a two-step cascade, we expect |0 that the 
I^rp spectrum would be softer than that shown in Fig. 1, so that SM backgrounds to the 
signal could be problematic. For this reason, we focus on the signal where one (or both) 
of the charginos decay leptonically. This stop signal, which consists of events with 1 or 2 
isolated leptons together with jet(s) and (from the undetected Zi's and neutrino(s)), has 
the same event topology as the canonical signal for top quarks. It is, therefore, essential to 
devise strategies to separate the stop signal events from top background, keeping in mind 
that cuts designed to optimize the top quark signal may not be suitable for the detection of 
ti squarks [|l^. 

A. The Single Lepton Signal. 

At the Tevatron, single lepton events from the cascade decays of gluinos and squarks are 
expected to be swamped by backgrounds from high px W production, with the isolated 
lepton coming from the decay of the W. What is different about stop events is that each 
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event contains two hard, central 6-quark jets which may be tagged using a sihcon microvertex 



detector. In our computation, we assume [19| 40% of the events have the vertex within the 



SVX barreL For these events, we take the efficiency for detecting i?-hadrons with > 15 
GeV and |?7b| < 1 to be 30%. In our analysis of the li signal, we have required the following: 

• (z) one isolated e (fi) with px > 10 GeV (5 GeV), 

• (a) two or more jets with at least one of the jets in the central region \rij\ < 2, 

• (m) I^j. > 25 GeV, and 

• (iv) at least one tagged 5-hadron. 

The main SM backgrounds to the single lepton signal come from W —>■ iu + jet production, 
and from tt production and are shown in Table I for = 140 GeV and rrit = 170 GeV along 
with the signal with m^-^ = 100 GeV, = 70 GeV and = 30 GeV. We emphasize 
that 6-tagging capability as assumed above is absolutely essential for the identification of 
the signal above the W background. ISAJET generates Wbb events from initial state gluon 
radiation followed by gluon splitting, using the usual coUinear QCD matrix elements of the 
parton shower model. We have checked that without the 6-tagging requirement (iv), the 
signal to W background ratio would be degraded by a factor > 250. 

We see from Table I that the stop signal is considerably smaller than the total SM back- 
ground. We have studied distributions that may serve to enhance the signal-to-background 
ratio. The multiplicity distribution of jets in the stop signal (solid) and in the top back- 
ground for nit = 140 GeV (dotted) and nit = 170 GeV (dashed-dotted), after cuts (i-iv) 
is shown in Fig. 5 a. We have illustrated the signal for m^^ = 100 GeV, = 70 GeV, 
and m^^ = 30 GeV. As expected, the top background has a significantly larger jet multi- 
plicity than the signal. We have checked that this is true even if ti is as heavy as 130 GeV. 
To enhance the stop signal to top background ratio without much loss of signal, we have 
required 

• (v) rijet < 4. 

In this choice we have been guided by the fact that the number of signal events in the 
data sample of the current Tevatron run, which is expected to accumulate an integrated 
luminosity of 100 pb~^, is not very large. It is clear, however, that if a significantly larger 
data sample is available (as should be possible after the main injector upgrade), choosing 
^jet < 3 (recall we always require at least two jets in the signal) optimizes the stop signal 
to top background ratio. 

The signal-to-background ratio can be further improved by recognizing that in both the 
top and W backgrounds the lepton and much of the come from the decay of a single 
W^-boson. We thus expect that the transverse mass distribution da/dmT{i, I^rp) for the 
background exhibits a characteristic Jacobian peak while the signal should show no such 
feature. Toward this end, we show this distribution, after cuts (i-v) in Fig. 56 for the same 
signal case as in Fig. 5a (solid), the W iu background (dashes) and the ti background 
for rut = 140 GeV (dotted). The transverse mass distribution will thus exhibit substantial 
distortion for low mxii, Ifj) if a signal is present. We see that by further requiring. 
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• (m) mTie,:^T) < 45 GeV, 

we are able to eliminate a large fraction of the background with a relatively modest loss of 
signal. The background, including cuts (v) and (vi), is shown in the last column in Table I. 
We see that a signal-to-background ratio of unity is possible for m^-^ < 100 GeV, even for 
the unfavourable = 140 GeV case. Finally, we examined whether tails from direct bb 
production, which has an enormous cross section at the Tevatron, could possibly fake the 
stop signal. We simulated 3M bb events, and found that only one event passes our cuts. 
Applying the 6-tagging efficiency, we estimate that this background to be ~ 24 fb, with 
large errors. 

In Fig. 6, we show contours of constant cross section for the single lepton signal in the 
vs. plane for n = —rrig and tan (3 = 2. The corresponding Zi mass is also shown on 



the right. We have cut off the contours for m^^ < 20 GeV since this region is excluded ||20 



within the MSSM framework. We see that unlike Fig. 3, where the cross section remains 
substantial even for Zi masses close to m^^, the single lepton signal rapidly decreases when 
™ii approaches rn^^. This is because the 6-jet becomes very soft and so is not tagged. We 
also see that ^ 20 signal events (compared to a background of 6-11 events, depending on nit) 
may be expected in a Tevatron data sample of 100 pb~^ for m^-^ < 90 GeV if fTi^^ ^ 60 GeV. 
This corresponds to about a 6a effect in the data sample expected to be accumulated by 
the end of the current Tevatron run, provided that sufficient b-tagging efficiency is achieved. 
In the longer term, the Tevatron may be able to probe a stop mass in excess of 100 GeV in 
this channel, particularly after the main injector becomes operational. 



B. The Dilepton Signal. 

Finally, we consider the dilepton signal arising from titi events where both the charginos 
decay leptonically. For this signal, we require: 

• (i) at /east two unlike sign, isolated leptons (e or /i), with pT^e) > 8 GeV andpr(/^) > 5 
GeV, 

• (zz) 20° < A(f){i+J'-) < 160° (to reduce Drell-Yan background), 

• (m) at least one jet with \T]j\ < 2, 

• (iv) > 25 GeV. 

The SM backgrounds to this signal, which dominantly come from ti production and W 
pair production, are shown in Table II after the cuts (i-iv). Also shown is the bound on 
this cross section we get from a simulation of lOM bb events, none of which pass our cuts. 
For comparison, we also show the signal for the same case as in Table I. We see that the 
signal is smaller than the SM backgrounds even for the heavier top case. Unlike the case 
for the single lepton signal, there is no obvious kinematic variable like the transverse mass 
that can be used to enhance the signal relative to the background. We note, however, that 
the leptons from the three-body chargino decays are expected to be softer than background, 
which has two-body decays of real W bosons. 
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With this in mind, we examined several distributions (to be discussed shortly) that may 
help to distinguish the top squark signal from the SM backgrounds. We found that the 
distribution of the scalar sum of the transverse momenta of the two leptons plus the I^j, in 
each event, which we refer to as "bigness", or S, 

B=\pr{t)\ + \pT{n\ + WTl (8) 

provides optimal distinction between the signal and the (dominant) top quark background. 
This distribution is shown in Fig. 7a for the signal from a 100 GeV stop decaying into a 
chargino (solid) and for the ti background for irtt = 140 GeV. As anticipated, the distribution 
is considerably harder for the background than for the signal; requiring 

• (v) B < 100 GeV, 

in addition to cuts {i-iv), significantly improves the signal to background ratio. The SM 
backgrounds, including this cut, are shown in the last column of Table II. The difference in 
the signal and background "bigness" distributions is illustrated in an alternative fashion in 
Fig. 7b, where we have shown on the horizontal (vertical) axis the signal (top background) 
cross section integrated up to the value of B marked on the curve in this figure. In other 
words, the signal (background) cross section is 80 fb (15 fb) if i? < 80 GeV. It is clear 
that the best distributions for distinguishing the signal from background are those where 
the slope of the corresponding curve remains small for a large distance after which it turns 
upward. Furthermore, since we can read off the value of B for any point along the curve. 
Fig. 7b enables us to easily find an optimal value for the cut. We see that B < 100 GeV 
that we obtained from the histogram in Fig. 7a is indeed a good choice for the cut. 

In Fig. 8, we have shown the analogous curves for other combinations of the lepton 
momenta and I^rp that may serve to distinguish the signal from background for the same 
cases as in Fig. 7. We show (a) the bigness distribution (again, for ease of comparison), 
(6) the I^j, distribution, (c) the distribution of the scalar sum, |pt(^^)| + \pt{^'~)\, and (d) 
the pt distribution of the hard lepton. From the shapes of these curves, it should be clear 
that the "bigness" distribution is indeed the best one to distinguish between stop and top 
production. 

We have illustrated the px distributions of the two leptons in the dilepton signal as well 
as the corresponding I^rp distributions, after cuts {i-v), in Fig. 9 for (a) m^^ = 100 GeV, 
m^^ = 90 GeV and m^^ = 42 GeV, (b) m^^ = 100 GeV, m^^ = 60 GeV and m^^ = 23 
GeV, and (c) m^^ = 70 GeV, = 60 GeV and m^^ = 23 GeV. The fast lepton and 
distributions in cases (a) and (6) are considerably harder than in case (c) where the stop 
is rather light and rn^^ — m^^ is small. Also, for a given value of m^^, these distributions 
become softer with decreasing difference between the chargino and the LSP mass. The slow 
lepton Pt distributions are, however, backed up against the cut. These distributions also 
enable the reader to assess the effect on the signal if the cut on the lepton pt is further 
relaxed. Clearly, capability for detection of leptons near the edge of the Pt(^) cut (z) in 
events triggered by jets and/or I^j, is required in order not to lose too much of the signal. 

In Fig. 10, we show the cross section contours for the dilepton signal after all the cuts for 
rrig = —fi, and tan/5 = 2 in the m^^ — plane. The region below the horizontal dashed 
line is experimentally excluded ||2^ within the MSSM framework, while in the region above 
the sloping dashed line, the stop decays via the flavour- changing loop decay mode analysed 
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in Sec. III. We see that while the signal is rather small, it exceeds the SM backgrounds for 
rrit^ < 110 GeV (130 GeV) for = 140 GeV (170 GeV) even if Zi is heavy. In a data 
sample with an integrated luminosity of 100 pb~^ that should be accumulated by the end of 
the current Tevatron run, we may expect > 10 signal events compared to a background of 3-6 
events, depending on the top quark mass, if m^^ < 100 GeV. Thus the Tevatron experiments 
should soon be able to probe top squark masses up to just under 100 GeV (larger, if the 
top is discovered and found to be heavy) even in the dilepton channel. This should increase 
to 120-130 GeV after the main injector becomes operational and the top turns out to be 
heavy. We remark that unlike as for the single lepton signal, the detection of the signal in 
the dilepton channel does not rely on the capability for tagging displaced vertices. 

We remind the reader that in all our computations up to this point, we have chosen 
parameters so that the leptonic branching fraction for chargino decays is 11% per lepton 
family. As mentioned above, this may be considerably enhanced |2l|] if sleptons are lighter 
than squarks and the Zi is mainly a U(l) gaugino. For instance, for rriq = rrig = 250 GeV, 
/i = —400 GeV with tan/3 = 20 (which yields = 70 GeV, m^^ = 35 GeV and slepton 
masses around 125 GeV), this branching fraction increases to about 30%. For this case, 
we find the 6-tagged single-lepton (dilepton) cross section after all cuts of 220 fb (780 fb) 
for a stop mass of 100 GeV to be compared with 140 fb (130 fb) in Table I. Our point 
here is to illustrate that the leptonic cross sections can indeed be significantly enhanced for 
reasonable choices of SUSY parameters, and that these enhancements are in keeping with 
expectations from the increase in the leptonic branching fraction of the chargino. In such 
scenarios, it may be possible to probe stop masses significantly beyond 100 GeV by the end 
of Run I B of the Tevatron. To illustrate, we have checked that for m^-^ = 130 GeV, the 
choice rrig = rriq = 315 GeV, /i = —500 GeV and tan/3 = 20, yields cross sections of 63 fb 
and 210 fb in the single-lepton and dilepton channels, respectively. It is amusing to see that 
the dilepton signal can be considerably larger than the signal in the single-lepton channel. 



V. DISCUSSION AND CONCLUSIONS 

We have seen that, in the framework of supergravity models, the lighter of the two 
top squarks (ti) may be very light even if the other squarks and gluinos have masses of 
several hundred GeV. Current Tevatron experimental analyses on squark masses are 
not applicable to a single light top squark, so that currently the best limit is still m^-^ > 45 
GeV, based on non-observation of a signal at LEP. Hence, one of most favorable routes 
towards the discovery of supersymmetry at the Fermilab Tevatron collider may be to search 
for the light top squark. We have shown in this paper that Tevatron experiments ought to 
be able to explore top squark masses of ~ 100 GeV, given an integrated luminosity of ~ 100 
pb~^. Thus, Tevatron collider experiments ought to be able to either discover supersymmetry 
via the top squark, or place an important new limit on its mass. Such a limit can serve to 
constrain the GUT scale parameter space of supergravity models, where the ti mass is driven 
to small values by a large top quark Yukawa coupling. In fact, model builders already have 
to take care to ensure that m?^ is not driven negative, leading to breaking of electromagnetic 
and colour gauge symmetries |]l|. 

If mi_^ < nib + '^fyj) then ii cZi, so collider experiments ought to search for pp 

titi —>■ ccZiZi, leading to a multi-jet plus I^j. signature. Even after suitable cuts, significant 
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backgrounds from W and Z production remain. The signal cross section, which depends 
only on the ti and Zi masses, is conveniently summarized by cross section contours in the 
— m^^ plane in Figs. 3 and 4. Assuming that the backgrounds can indeed be reliably 
suppressed/subtracted out from the data sample (via r-jet veto, Z+jets total normalization, 
etc.), we see that ii masses of ~ 80 — 100 GeV ought to be accessible to Tevatron experiments 
by the end of the current run. 

If m^-^ > rrib + m^^, then ti bWi, so collider experiments ought to search for pp 

iiii — » bbWiWi, where Wi iuZi or qq'Zi. Just as for the top quark search, signals 
appear most promising in the single and dilepton final states. For single lepton final states, 
5-tagging capability is essential to veto the large M/^+multi-jet background. The results of 
our computation of the single-lepton and the dilepton signal are summarized in Fig. 6 and 
Fig. 10, respectively. Given sufficient 5-tagging capability, ii masses to 80-100 GeV ought to 
be explorable at the Tevatron via the single lepton channel. For the dilepton signal, major 
backgrounds come from ti and WW production. Requiring "soft" dilepton events, [e.g., 
requiring B = |pt(^i)| + \pt{^2) \ + -^t < 100 GeV), allows sufficient rejection of background 
that, as before, stop masses of ~ 100 GeV ought to be probed at the Tevatron by the end 
of the current run. 

The analysis in the present paper is a conservative one, in that an additional source 
of stop squark production, t ^i^i, has been neglected. This decay mode, which can 
have branching fractions of up to ~ 30%, would contribute an additional signal, as well as 



diminish the SM top quark background |]T5[. It would be interesting to investigate whether 
the observation of a top quark signal at a rate consistent with SM expectations could be 
used to independently constrain m^-^ at a top factory. We have also been conservative in 
that for the leptonic signals from stop production shown in Figs. 6 and 10, we have chosen 
parameters so that the branching fractions for leptonic decays of charginos are essentially 
the same as those of the W boson; as discussed in the text, the single lepton (dilepton) 
signals may be enhanced by a factor of ~1.5 (~6) for certain ranges of SUSY parameters. 

In summary, the search for stop squarks adds a new avenue to the search for supersym- 
metry at the Tevatron collider. Light top squarks have a significant degree of theoretical 
motivation. We have attempted to show in this paper that our experimental colleagues 
ought to be able either to discover top squarks at the Tevatron, or to significantly extend 
the limit on their mass. 
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TABLES 



TABLE I. Cross sections in fb at the Tevatron after cuts described in the text for SM back- 
ground to the S-tagged single lepton signal. Also shown is the signal cross section for a represen- 
tative case near the edge of the Tevatron reach, where we have taken m^^ = 100 GeV, m^^^ = 70 
GeV and m-x = 30 GeV. Results have summed over e and /li. 



process Cuts (i-iv) cuts{i-vi) 

W ^iu 440 30 

tt(140) 450 81 

tt(170) 170 32 

bb_ 24 24 

tifi(lOO) 240 140 



TABLE II. Cross sections in fb at the Tevatron after cuts described in the text for SM back- 
grounds to the dilepton signal from top squark production. Also shown is the signal cross section 
for the same case as in Table I. Results have summed over e and 



process Cuts (i-iv) cuts{i-v) 

W+W- 50 10 

tt{UO) 360 54 

tt(170) 135 14 

bb_ < 48 

fifi(lOO) 180 130 
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FIGURES 



FIG. 1. ^j, distributions, before any cuts, from tiii production in pp collisions at -y/s = 1.8 TeV 
for two choices of stop and Zi masses. Here, the stop is assumed to decay via t\ cZi. Also 
shown are SM background distributions from Z ^ vv production (dotted), W ^ tv production 
(dashed) and W ^ [iv production (dot-dashed). 

FIG. 2. The scatter plot for events in the A(l){I^rp, j^ear) - PT{jfast) plane after cuts {i-iv) 
in Sec. Ill of the text for the ti^i signal at the Tevatron for (a) = 85 GeV, m-^^ = 20 GeV, 
{h) mi_^ = 125 GeV, m^^ = 40 GeV, and dominant SM backgrounds from (c) W ^ tu production, 
and (d) Z ^ vv production. 

FIG. 3. Cross section (in ph) contours for the distributions, after cuts [i-v) in Sec. Ill of 
the text, from t-^x production in pyp collisions at s/s = 1.8 TeV, assuming that the stops decay via 

ti — > cZi. Also shown are dominant backgrounds from SM processes. The W ^ t background 
will be smaller if it is possible to veto r jets at some level. Also, it may be possible to subtract the 
Z ^ vv background as discussed in the text. 

FIG. 4. The same as Fig. 3 except that, in this figure, we require a non-isolated muon with 
PrifJ-) > 3 GeV to tag the c-jet. Cuts are as described in Sec. Ill of the text. 

FIG. 5. (a) The jet multiplicity distribution in S-tagged single lepton events from stop pair 
production at a pp collider at ^/s = 1.8 TcV, assuming that stop decays via ti bWi after cuts 
(i-iv) in Sec. IV. A of the text. We illustrate the distribution from the signal (solid) for m^^ = 100 
GeV, Tn^^ = 70 GeV and = 30 GeV, and from top quark backgrounds for mt = 140 GeV 
(dotted) and = 170 GeV (dashed-dotted) , and (6) the transverse mass 771^(^,-^7^) distribution 
in these events after cuts (i-v) in Sec. IV.A of the text. We illustrate the distribution from the 
signal (solid) for the same SUSY masses as in (a), and from the dominant SM backgrounds from 
W ^ £u production (dashed) and ti production for mt = 140 GeV (dotted). The branching 
fraction for lUi evZ\ decay is 11%. If the top quark has a mass of 170 GeV, its contribution to 
the background is significantly smaller. 

FIG. 6. Cross section (in jh) contours in the — plane for S-tagged single lepton 

events from stop pair production at a pp collider at ^fs = 1.8 TeV, assuming that stop decays via 
ti bWi after cuts (i-vi) in Sec. IV.A of the text. The branching fraction for Wi evZ\ decay 
is 11%. Also shown are the main SM backgrounds. The mass of Z\ is illustrated on the right. 
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FIG. 7. (a) Distributions of the "bigness" variable B defined in Sec. IV.B of the text for 
dilepton events from stop pair production at the Tevatron (where the stop decays via ti bWi) 
and from the background from production of 140 GeV ti pairs after cuts {i-iv) in Sec. IV.B of the 
text. The sparticle masses and chargino branching fractions are as in Fig. 5. (&) The cross section 
for stop dilepton signal (top background) on the horizontal (vertical) axis, integrated up to the 
value of B on the curve after the same cuts as in (a). This plot is described in more detail in Sec. 
IV.B. 

FIG. 8. Cross sections with (a) the bigness, (6) I^j,, (c) |pt(-^i)| + \pt{'^2)\, and (d) pxi^fast) 
variables integrated up to the values shown on the respective curves, for stop dilepton events 
(horizontal axis) and tt background, for nit = 140 GeV, after cuts (i-iv) of Sec. IV.B. 

FIG. 9. Distributions of the transverse momentum of the two leptons in dilepton events from 
stop pair production at the Tevatron along with the distribution after cuts (i-v) in Sec. IV.B 
for three illustrative cases of ii, W\ and Zi masses mentioned in the text. 

FIG. 10. Cross section (in fh) contours for the dilepton cross section from the production of 
stop pairs at a pp collider with yfs = 1.8 TeV, including cuts {i-v) of Sec. IV.B. The branching 
fraction for W\ euZi decay is 11%. Also shown are the main SM backgrounds. 
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